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ABSTRACT 

The objective of this paper is to determine the level of obscured star formation activity and dust 
attenuation in a sample of gamma-ray burst (GRB) hosts; and to test the hypothesis that GRB hosts 
have properties consistent with those of the general star-forming galaxy populations. We present a radio 
continuum survey of all z < 1 GRB hosts in The Optically Unbiased GRB Host (TOUGH) sample 
supplemented with radio data for all (mostly pre-Swift) GRB-SN hosts discovered before 2006 October. 
We present new radio data for 22 objects and have obtained a detection for three of them (GRB 980425, 
021211, 031203; none in the TOUGH sample), increasing the number of radio-detected GRB hosts from 
two to five. The star formation rate (SFR) for the GRB 021211 host of ~ 825 M Q yr~\ the highest ever 
reported for a GRB host, places it in the category of ultraluminous infrared galaxies. We found that at 
least ~ 63% of GRB hosts have SFR < 100 M yr" 1 and at most - 8% can have SFR > 500 M Q yr" 1 . 
For the undetected hosts the mean radio flux (< 35/iJy 3a) corresponds to an average SFR < 15 M 
yr _1 . Moreover, > 88% of the z < 1 GRB hosts have ultraviolet dust attenuation A\jy < 6.7 mag (visual 
attenuation Ay < 3 mag). Hence we did not find evidence for large dust obscuration in a majority of 
GRB hosts. Finally, we found that the distributions of SFRs and A\jy of GRB hosts are consistent 
with those of Lyman break galaxies, Ha emitters at similar redshifts and of galaxies from cosmological 
simulations. The similarity of the GRB population with other star-forming galaxies is consistent with 
the hypothesis that GRBs, a least at z < 1, trace a large fraction of all star formation, and are therefore 
less biased indicators than once thought. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: ISM — galaxies: star formation 
— gamma-ray burst: general — radio continuum: galaxies 



1. INTRODUCTION 

Long gamma-ray bursts (GRBs) mark the endpoint of 
the lives of very massive stars (e.g. Hjorth et al. 2003; 
Stanek et al. 2003) and due to the short life-times of such 
stars, they are believed to be excellent tracers of ongoing 



star formation in distant galaxies (Jakobsson et al. 2005; 
Yiiksel et al. 2008; Kistler et al. 2009; Butler et al. 2010; 
Elliott et al. 2012; Robertson & Ellis 2012). However, be- 
fore GRBs can be quantitatively used to trace the star 
formation history of the Universe, the properties of their 
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host galaxies and the biases of the GRB samples must be 
understood. 

From optical/near- infrared studies we know that GRB 
hosts are often faint dwarf galaxies (Le Floc'h et al. 2003; 
Christensen et al. 2004; Savaglio et al. 2009; Castro Ceron 
et al. 2010; Levesque et al. 2010b; Svensson et al. 2010). 
However, some host galaxies of (often optically-obscured) 
GRBs are massive (M* > 1O 1O ' 5 M ) and/or belong to 
the category of luminous infrared galaxies (LIRGs; Lir > 
lO n L , or star formation rate SFR > 17.2 M Q yr 1 us- 
ing the conversion of Kcnnicutt 1998), e.g. GRB 980613 
(Castro Ceron ct al. 2006, 2010), 020127 (Bcrgcr ct al. 
2007), 020819B (Savaglio et al. 2009; Kiipcii Yoldas, ct al. 
2010), 051022 (Castro-Tirado ct al. 2007; Savaglio et al. 
2009), 070306 (Jaunsen et al. 2008; Kriihler et al. 2011), 
080207 (Hunt et al. 2011; Svensson et al. 2012), 080325 
(Hashimoto et al. 2010), and 080605 (Kriihler et al. 2012a). 
Hence, the diversity of the GRB host sample is not yet fully 
described. 

Moreover, short- wavelength emission does not give us a 
complete picture of GRB hosts, as it misses star forma- 
tion that is heavily obscured by dust. Unfortunately long- 
wavelength emission has been detected only in a handful 
of GRB hosts (Berger et al. 2001, 2003a; Frail et al. 2002; 
Tanvir et al. 2004; Castro Ceron et al. 2006; Le Floc'h et al. 
2006, 2012; Priddey et al. 2006; Michalowski et al. 2009; 
Stanway et al. 2010; Hunt et al. 2011; Watson et al. 2011; 
Hatsukade et al. 2012; Svensson et al. 2012; Walter et al. 
2012; see the compilation of submillimeter observations in 
de Ugarte Postigo et al. 2012), which hampers our ability 
to study GRB hosts in the context of galaxy evolution, as 
a significant fraction of star formation in the Universe is 
believed to be obscured by dust (e.g. Lc Floc'h ct al. 2005; 
Chapman et al. 2005; Wall et al. 2008; Michalowski ct al. 
2010a). 

Here we attempt to improve this situation through an 
extensive multi-facility radio survey of GRB hosts limited 
to z < 1 to obtain meaningful SFR limits, drawn from 
The Optically Unbiased GRB Host (TOUGH; Hjorth et al. 
2012) survey, which allows for unbiased statistical analysis. 

Observations at radio wavelengths provide an unob- 
scured (unaffected by dust) view on star-forming galaxies 
by tracking directly the recent (< 100 Myr) star formation 
activity through synchrotron radiation emitted by rela- 
tivistic electrons accelerated by supernova (SN) remnants 
(Condon 1992). Moreover, even though the radio emission 
accounts for only a fraction of the bolomctric luminosity 
of a galaxy, it is well correlated with the infrared emis- 
sion, a good tracer of both the SFR and the dust mass in 
a galaxy. Finally, the timescale probed by radio emission 
(< 100 Myr) is much longer than the lifetime of a GRB 
progenitor (~ 5-8 Myr; Sollerman et al. 2005; Hammer 
et al. 2006; Ostlin et al. 2008; Thonc et al. 2008), so the 
radio emission probes the average star formation state of 
a galaxy, unlike a GRB rate, which measures the almost 
instantaneous star formation activity. 

The objective of this paper is to (1) determine the level 
of obscured star formation activity and dust attenuation in 
a representative sample of z < 1 GRB hosts, and (2) test 
the hypothesis that GRB hosts are consistent with the gen- 
eral star-forming galaxy populations at similar redshifts. 

19 http : //www . dark-cosmology . dk/TOUGH 



We use a cosmological model with Hq = 70 km s _1 
Mpc -1 , Ha = 0.7, and f2 m = 0.3 and assume the Salpctcr 
(1955) initial mass function to which all estimates from 
the literature were converted to, if necessary. 

2. SAMPLE 

Our target sample is composed of two subsets. The 
main subset is drawn from the TOUGH sample based on 
the Swift satellite and a Very Large Telescope (VLT) Large 
Programme. The survey design, the selection criteria, and 
the summary of the host properties (including redshifts) 
are presented in Hjorth et al. (2012), the photometry and 
the host properties are analyzed in Malcsani et al. (2012), 
the redshifts are presented in Jakobsson et al. (2012) and 
Kriihler et al. (2012b), and the Lya properties are dis- 
cussed in Milvang- Jensen et al. (2012). The reduced data 
will be available from the TOUGH Web site 19 . The sam- 
ple includes all long GRBs that exploded between 2005 
Marchl and 2007 August 10, observable from the south- 
ern hemisphere (—70° < 6 < +27°), with low Galactic 
foreground extinction (Ay < 0.5 mag) and no bright star 
nearby, for which X-ray observations are available < 12 
hr after the burst (with < 2" error circle radius) to al- 
low the determination of accurate positions. Therefore, 
this X-ray-selected sample is constructed in a way that 
it is not biased against dusty systems and the selection 
docs not depend on the host luminosity. We note that 
the availability of redshift does depend on the host lu- 
minosity, but redshifts were measured for ~ 77% (53/69) 
of Swift/VLT TOUGH GRBs (Hjorth et al. 2012; Jakob- 
sson ct al. 2012). Moreover, half of the TOUGH GRB 
redshifts were obtained from optical observations of after- 
glows, so the redshift recovery fraction is dependent on 
the host brightness only for the other half. Finally, fainter 
hosts (for which redshifts could not be measured) are less 
likely to be at z < 1. Indeed 75% (12/16) of TOUGH 
GRBs with unknown redshifts are fainter than R = 25 
mag, whereas the same is true only for ~ 17% (2/12) of 
GRBs confirmed to be at z < 1. We therefore conclude 
that the z < 1 TOUGH sample has a completeness nearing 
100%, and, in any case, larger than ~ 77%. 

We restricted the TOUGH sample to z < 1 to obtain 
meaningful radio constraints on SFRs. The z < 1 TOUGH 
unbiased subset consists of 12 hosts of which all were ob- 
served within our program (see Table 1 for observation 
logs). GRB 060814 at z = 1.92, 050915A at z = 2.527, 
and 070808 with currently no redshift measurement are in 
our target sample, because they were initially believed to 
be at z < 1 (e.g. Thone et al. 2007). 

The second subset includes (mostly pic-Swift) GRBs 
that were spectroscopically or photometrically confirmed 
to be associated with SNe before 2006 October, namely 
the sample of Fcrrcro et al. (2006) plus GRB 980425/SN 
1998bw (Galama et al. 1998) and GRB 040924 (Sodcr- 
bcrg et al. 2006). We targeted GRB-SN hosts, because 
their progenitors are securely established to be connected 
with recent star formation (see Hjorth & Bloom 2011 for 
a recent review of the GRB-SN connection). Since the de- 
tection of an SN component in a fading GRB afterglow is 
difficult at high redshifts, this selection imposes a practi- 
cal limit of z < 1. In total 15 hosts were selected (with 



Radio survey of GRB hosts 



3 



Table 1 
Radio Observation Logs 



GRB 


Array 


Observation Dates 


t" 


Frequency 


rms 


Synth. Beam Size 


Calibrators 








Or) 


(GHz) 


(MJy) 


(") 












GRB-SN 


subset 






980425 c d 


ATCA 


2007 Aug 18 


9.00 


4.8, 8.64 


46. 27 


76 X 38, 37 x 21 


PKS B1934-638 


991208 


WSRT 


2007 Aug 2-3 


11.97 


1.43 


47 


14.5 X 10.5 


3C286, 3C48 


020903 


GMRT 


2008 Jan 18-19, Mar 1 


16.99 


1.43 


41 


4.0 x 2.0 


3C48, 2243-257 


021211 


VLA 


201)7 Jul 14 


5.45 


1.43 


31 


1.7 x 1.5 


Iao 0542+498, PMJN J0808+0514 


031203 d " 


ATCA 


2008 Jan 26 


6.97 


1.39, 2.37 


46, 37 


8.5 x 3.4, 6.3 x 2.3 


PKS B1934-638, PKS B0826-373 


U41UUO 


CjMKl 


2007 Aug 7—8 


9.61 


1.43 


181 


7.3 x 2.0 


3C48, J32 0026+34 








TOUGH z < 1 unbiased 


su bsct 




n t \ 1 1 c a 
UOU1J.DA 




/UUo /\pr z/— zo 


11.97 


1.43 


75 


oi.u x y.o 


Oli4fi, oLjZOU 




VV ortl 


ZUU i Aug lo — 


11.97 


1.43 


52 


oo. / X 141.5 




050824 


WSRT 


2007 Doc 26-27 


11.97 


1.43 


100 


39.9 X 15.5 


3C147, 3C286 


051016B 


WSRT 


2007 Dec 28-29 


11.97 


1.43 


47 


59.7 X 14.1 


3C147, 3C286 


05iH7B rf 


ATCA 


2009 Aug 12 


8.29 


5.5. 9.0 


12, 19 


6.4 x 1.7 3.9 x 1.1 


PKS B1934-638 PKS B0607-157 


06021 8 


WSRT 


2007 Aiip- 1 6—1 7 


11.96 


1.43 


117 


51.9 x 15.2 


3C!48 3C286 


06061 4 rf 


ATCA 


2009 A hp- 9-10 


8.84 


5.5, 9.0 


11, 14 


3.1 x 1.9. 1.7 x 1.0 


PKS B1 934-638 


060729 


ATCA 


2008 Jan 26,28 


11.36 


1.39 


35 


7.4 x 6.4 


PKS B1934-638, PKS B0515-674 


060912A / 


GMRT 


2009 Jun 1-2 


10.00 


1.43 






3C48 


061021 


ATCA 


2008 Apr 18 


7.90 


1.39 


36 


20.0 x 4.8 


PKS B1934-638, PKS B0919-260 


oemoA^ 


WSRT 


2007 Dec 29 


12.00 


1.43 






3C147, 3C286 


070318 


ATCA 


2008 Apr 19 


9.74 


1.39 


47 


7.2 x 4.2 


PKS B1934-638, PKS B0405-385 










Other 


hosts 






050915A 


ATCA 


2008 Jan 25,27 


15.62 


1.39 


29 


18.3 x 5.5 


PKS B1934-638, PKS B0451-282 


d 


ATCA 


2011 Dec 19 


9.78 


5.5, 9.0 


12, 15 


5.9 x 2.1, 3.7 x 1.3 


PKS B1934-638, PKS B0537-286 


060505 d 


ATCA 


2009 Aug 10-11 


8.48 


5.5, 9.0 


17. 14 


5.2 x 1.7, 3.4 x 1.1 


PKS B1934-638, PKS B2155-152 




GMRT 


2008 Jan 20-21 


5.89 


1.43 


58 


3.6 x 2.3 


3C48, 2243-257 


060814 


WSRT 


2007 Dec 30 


11.96 


1.43 


7M 


42.7 X 14.9 


3C48, 3C286 


070808 


GMRT 


2009 Jun 2-3 


6.87 


1.43 


68 


3.6 x 2.2 


3C48, 0022+002 



Note. — The horizontal lines divide the GRB-SN and the z < 1 TOUGH subsets (see Section 2) and the hosts which do not belong to any 
of these subsets. GRBs 050525A and 060218 belong to the first two subsets. 
a On-source integration time. 

6 The first (second) object was used as a primary (secondary) calibrator. For WSRT the indicated objects were used as primary calibrators at 
the beginning and the end of the run. 
c Data published in Michaiowski et al. (2009). 

d This object was observed simultaneously at two frequencies, see Table 2. 
e Data published in Watson et al. (2011). 

f Poor quality (interference and system malfunctions) of the data impedes the flux density measurement. 



an overlap of two hosts with the TOUGH subset) of which 
eight were observed within our program and for the re- 
maining seven the deep radio upper limits from the liter- 
ature were adopted (see Table 2). 

In summary, our sample consists of 30 GRB hosts and 
we provide new radio observations for 22 of them. 

3. DATA 

The radio data were collected using the Australian Tele- 
scope Compact Array (ATCA; proposals C1651, C1741, 
CX228) in the 6 km configuration (H168 for GRB 980425), 
the Giant Metrewave Radio Telescope (GMRT; propos- 
als 12MMc01, 13MMc01, 16_093), the Very Large Ar- 
ray (VLA; proposal AM902) in A configuration, and the 
Wcstcrbork Synthesis Radio Telescope (WSRT; proposals 
R07B004, R08A002) in maxi-short configuration. The log 
of observations is presented in Table 1. 

Data reduction and analysis were done using the 
MIRIAD (Sault & Killeen 2004) and AIPS 20 packages. 
Calibrated visibilities were Fourier transformed using "ro- 
bust" or "uniform" weighting depending on which gave a 
better result for a particular field. The resulting rms noise, 
beam sizes, and calibrators are listed in Table 1 and the 
radio contours overlaid on the optical images for detected 
hosts are presented in Figure 1. The data for the observa- 

20 http : //www . alps . nr ao . edu/ cook . html 



tions of GRB 060912A and 061110A were found to have 
been severely affected by radio frequency interference and 
system malfunctions. Therefore, we had to discard a sig- 
nificant fraction of the data and the remaining data were 
insufficient to create reasonable radio images of the fields. 

Flux densities were measured by fitting two-dimensional 
Gaussian functions to the region around the host and the 
errors were determined from the local rms on the images. 
The hosts of GRBs 980425 and 031203 slightly overlap 
with radio objects ~ 70" south (see Michaiowski et al. 
2009) and ~ 6" northwest, respectively, so their flux den- 
sities were estimated by simultaneous fitting two Gaussian 
functions with their centroids, sizes and orientations as free 
parameters. The lack of residuals left after the subtraction 
of these two Gaussians rules out a significant contamina- 
tion of the nearby objects to the measured flux densities 
of the hosts. 

4. RESULTS 

Our photometry measurements are presented in Table 2. 
Three (GRB 980425, 021211, 031203) out of twenty tar- 
geted hosts were detected (not counting upper limits from 
the literature) . Two out of these detections are in fact the 
first- and third-closest GRBs in our sample. None of the 
hosts in the TOUGH subset has been detected. Hence, 
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Table 2 

Radio Fluxes, Star formation Rates and Dust Attenuations of GRB Hosts 



GRB 


z 


Ref 


Flux Density 


Frequency 




Ref 


SFR radlo a 


SFRuv" 


Rcf c 










OJy) 


(GHz) 






(M yr" 1 ) 


(M e yr" 1 ) 




(mag) 


GR.B-SN subset 


970228 


0.695 


1 


< 69 


1.43 




30 


< 72 


0.60 


42 


<2.3 


980425 


0.0085 


2 


420 ± 50 


4.80 


+ 


31 


0.23 ± 0.02 


0.39 


42 


~ 








< 180 


8.64 


%, 


31 


< 0.17 


0.39 


42 


~ 


990712 


0.4337 


3; 4 


< 105 


1.39 




32 


< 36 


1.28 


43 


<1.6 








< 36 


5.50 




33 


< 35 


1.28 


43 


<1.6 








< 129 


9.00 




33 


< 180 


1.28 


43 


<2.4 


991208 


0.7063 


5 


< 32 


1.43 




i 


< 35 


0.83 


43 


<1.8 


000911 


1.058 


6 


< 57 


8.46 




34 


< 608 


1.40 


42 


<3.0 


010921 


0.451 


7 


< 83 


1.43 




30 


< 32 


1.60 


42 


<1.5 


011121 


0.36 


8 


< 120 


4.80 




30 


< 68 


1.83 


44 


<1.8 


020405 


0.691 


9 


< 42 


8.46 




35 


< 165 


3.70 


42 


<1.9 


020903 


0.251 


10 


< 53 


1.43 




X 


< 5.39 


0.42 


44 


<1.3 


021211 


1.006 


11 


330 ± 31 


1.43 




X 


825 ± 77 


0.72 


42 


3.4 








< 34 


2.10 




36 


< 114 


0.72 


42 


<2.5 








< 45 


8.46 




37 


< 427 


0.72 


42 


<3.1 


030329 


0.168 


12 


< 420 


1.40 




38 


< 17 


0.14 


42 


<2.4 


031203 


0.105 


13 


254 ± 46 


1.39 


t: 


39 


3.83 ± 0.69 


4.30 


12 


~ 








191 ± 37 


2.37 


t: 


30 


4.29 ± 0.83 


4.30 


42 


~ 








216 ± 50 


5.50 




33 


9.13 ± 2.11 


4.30 


12 


0.4 








< 48 


9.00 




33 


< 3.09 


4.30 


42 


~ 


040924 


0.859 


14 


< 63 


4.90 




40 


< 274 


0.66" 


X, 40 


<2.9 


041006 


0.716 


14 


< 45 


2.10 




36 


< 67 


0.47 


44 


<2.4 








< 348 


1.43 




X 


< 392 


0.47 


44 


<3.3 








< 123 


8.46 




41 


< 525 


0.47 


44 


<3.4 








TOUGH 


z < 1 unbiased 


subset 










050416A 


0.6528 


15; 16 


< 447 


1.43 




X 


< 405 


0.89 


44 


<3.0 


050525A 


0.606 


17 


< 228 


1.43 




X 


< 174 


0.64 


42 


<2.7 


050824 


0.828 


18; 19 


< 111 


1.43 




X 


< 175 


1.37 


45 


<2.4 


051016B 


0.9364 


20 


< 220 


1.43 




X 


< 465 


5.69 


X, 46 


<2.2 


051117B 


0.481 


21 


< 36 


5.50 




X 


< 44 


2.72 


t, 46 


<1.4 








< 57 


9.00 




X 


< 101 


2.72 


X, 46 


<1.8 


060218 


0.0334 


22 


< 447 


1.43 




X 


< 1.00 


0.05 


42 


<1.4 








< 117 


5.50 




33 


< 0.78 


0.05 


42 


<1.3 








< 42 


9.00 




33 


< 0.48 


0.05 


42 


<1.1 


060614 


0.125 


23 


< 33 


5.50 




X 


< 2.35 


0.02 


42 


<2.4 








< 42 


9.00 




X 


< 3.72 


0.02 


42 


<2.6 


060729 


0.54 


24 


< 105 


1.39 




X 


< 60 


0.13 


47 


<3.0 


061021 


0.3463 


25 


< 108 


1.39 




X 


< 22 


0.03 e 


t, 21 


<3.2 


070318 


0.836 


26 


< 141 


1.39 




X 


< 223 


1.44 


X, 48 


<2.5 


Other hosts 


050915A 


2.527 


21; 27 


< 59 


1.39 




X 


< 1204 


9.48 


X, 48 


<2.4 








< 44 


5.50 




X 


< 2521 


9.48 


X, 48 


<2.7 








< 37 


9.00 




X 


< 3032 


9.48 


X, 48 


<2.8 


060505 


0.0889 


28 


< 37 


5.50 




X 


< 1.50 


1.90 


42 


~ 








< 52 


9.00 




X 


< 2.53 


1.90 


42 


<0.1 








< 63 


1.43 




X 


< 1.05 


1.90 


42 


~ 


060814 


1.92 


21; 27; 29 


< 430 


1.43 




X 


< 4823 


31.20 


t, 48 


<2.5 


070808 


unknown 




< 156 


1.43 




X 











Note. — The horizontal lines divide the GRB-SN and the z < 1 TOUGH unbiased subsets (see Section 2) and the hosts which do not belong 
to any of these subsets. GRBs 050525A and 060218 belong to the first two subsets. For non-detected targets 3a limits are reported. 
"Assuming radio spectral index a = —0.75 and applying the calibration of Bell (2003). 
''From UV continuum unless noted otherwise. Not corrected for dust attenuation. 

c The symbol \ indicates that we derived the SFR from fluxes reported in the reference using the calibration of Kennicutt (1998) or Savaglio 
et al. (2009). 

d Visual extinction calculated from the ultraviolet extinction Atjv = 2.5 log(SFR radio /SFRTjv) assuming an SMC extinction curve, which gives 
A v = Auv/2.2. 
e From the [O II] line. 
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Fig. 1. — Radio contours (blue lines) overlaid on the optical images of the detected GRB hosts. The size of each imago depends on the host 
galaxy size and the resolution of the radio data: 3' for GRB 980425 (4.8 GHz), 10" for GRB 021211 (1.43 GHz), and 20" for GRB 031203 
(1.39 GHz). The red circles (with arbitrary sizes) mark the position of GRBs (optical positions for GRB 980425 and 021211 from Fynbo et al. 
2000; Fox et al. 2003, and X-ray position for GRB 031203 from Watson et al. 2004). The contours are 3, 4, 6, 8 and 10cr (see Table 1). North 
is up and east is left. The optical data arc from Sollerman et al. (2005, GRB 980425), Delia Valle et al. (2003, GRB 021211), and Mazzali 
et al. (2006, GRB 031203). 



this program (Michalowski et al. 2009, Watson et al. 2011, 
and this paper) increases the number of the radio-detected 
GRB hosts from two (GRB 980703, Berger et al. 2001; 
GRB 000418, Berger et al. 2003a) to five. Recently the 
host of GRB 031203 has also been detected at 5.5 GHz by 
Stanway et al. (2010). 

We assume that the entire flux is due to star formation 
and not active galactic nucleus (AGN) activity, which is 
a well-tested hypothesis for GRB hosts (see discussion in 
Michalowski et al. 2008; Watson et al. 2011). 

The SFRs derived from our radio data as well as from 
the ultraviolet (UV) data are presented in Table 2 and are 
shown as a function of redshift on Figure 2. The radio 
SFRs (SFR ra dio) were calculated from the empirical for- 
mula of Bell (2003, see Section 4.2 of Michalowski et al. 
2009 for discussion of its applicability to GRB hosts) as- 
suming a radio spectral index 21 a = —0.75 (Condon 1992; 
Ibar et al. 2010). This choice of spectral index has rela- 
tively small impact on derived SFRs, because our observed 
1.4 GHz data probe close to the rest-frame 1.4 GHz, at 
which the flux-SFR conversion is calibrated. Namely if we 
assumed a fiat index a = 0, then we would obtain SFRs 
~ 25-40% lower at z = 0.5-1. On the other hand, if we 
assumed a steeper value a = —1 (or —1.5), then we would 
obtain SFRs - 10-20% (~ 35-70%) higher at z = 0.5-1. 

The limit on the SFR of the GRB 980425 host based 
on 8.64 GHz data is not consistent with the value derived 
from the 4.80 GHz data, because, for consistency, a spec- 
tral slope of a = —0.75 was assumed, whereas in reality it 
is steeper (see Michalowski et al. 2009 and section 5.5). 

In order to assess amount of the dust attenuation in 
GRB hosts we compared their SFRs derived from the UV 
emission (SFRtjv) with SFR rac ii - In Table 2 we compiled 
the SFRtjv (mostly from 0.28 /jm continuum data) from 
the literature (Castro Ceron et al. 2010; Savaglio et al. 
2009; Christensen et al. 2004; Jakobsson et al. 2012; Ovald- 
sen et al. 2007; Svensson et al. 2010). The dc-reddened 

21 Defined as F v oc u a , i.e. a v v \ = log^^)/^^)]/ log(i^/n). 



SFRs given by Savaglio et al. (2009) were reddened based 
on their reported Ay. For the hosts of GRB 051016B, 
051117B, 060814 and 070318 we calculated the SFR UV 
from V-, and i?-band fluxes, respectively, reported by 
Ovaldsen et al. (2007) and Malesani et al. (2012), which 
correspond to the rest-frame UV emission at the redshifts 
of the hosts. For the hosts of GRB 040924 and 061021 
there are no UV continuum data available, so we calcu- 
lated SFR,q j T j from the flux reported by Wiersema et al. 

(2008) and Jakobsson et al. (2012), respectively, applying 
the conversion of Kewley et al. (2004, their equation (4)). 

We assume that SFR ra di reflects the total amount of 
star formation in GRB hosts. Hence, an approximate es- 
timate of the dust attenuation in the ultraviolet may be 
obtained by dividing the radio SFR and SFRtjv : 

A vv = 2.5 log S g ^^° ma S (!) 

The resulting attenuations are presented in Table 2. The 
uncertainties of SFR ra di and SFRtjv ar e of the order of 
a factor of two (Bell 2003; Kennicutt 1998), so the uncer- 
tainties of the Atjv estimates are of the order of a factor 
of 2^2 - 2.8 (~ 1.1 mag). 

5. DISCUSSION 

5.1. The ULIRG Nature of the Host of GRB 021211 

Our 1.43 GHz detection of the host of GRB 021211 
(~ lOer) corresponds to SFR of - 825M y r - 1 (Table 2), 
which places it in the category of ultraluminous infrared 
galaxies (ULIRGs; L m > 10 12 i Q , or SFR > 172M yr- 1 
using the conversion of Kennicutt 1998). This is the high- 
est SFR ever reported for a GRB host (compare with 
Berger et al. 2003a, 2001; Michalowski et al. 2008; Stanway 
et al. 2010; Watson et al. 2011). 

Because of this unusually high SFR we present the inves- 
tigation of the data quality for this object. We verified that 
the source is not due to an uncleaned bright source nearby. 
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Fig. 2. — Star formation rates (SFRs) as a function of redshift of GRB host galaxies. Squares and arrows denote SFRs derived from radio 
detections and 3<r upper limits, respectively. Circles denote lower limits on SFRs derived from the ultraviolet (UV) data. For a given GRB, 
the radio and UV SFRs are connected by a dotted line. GRBs are color-coded depending on whether they belong to the z < 1 TOUGH 
unbiased subset (red symbols), the GRB-SN subset (blue symbols), both (green symbols), or none (black symbols). The right j/-axis gives the 
corresponding infrared luminosity according to SFR(Mq yr _1 ) = 1.72 X 10 — 10 Lir(Lq) (Kennicutt 1998). The three low-redshift hosts (GRB 
980425, 031203, and 060505) are consistent with no dust attenuation as their SFR ra< ji are similar to SFRtjv- On the other hand, huge dust 
attenuation must be invoked to explain a very high SFR ra( ji of the host of GRB 021211. Crosses and plus symbols indicate the 24 (im and 
UV SFRs of Lyman break galaxies (LBGs; Basu-Zych et al. 2011) and Ho emitters (HAEs; mean values with standard deviations are shown; 
Sobral et al. 2009). GRB hosts are consistent with these populations. 



Moreover, the astrometry of our VLA map and the VLT 
image (Delia Valle et al. 2003) are consistent. Namely, for 
three radio sources that are detected in the optical image 22 
we measured the mean offsets between the radio and op- 
tical positions consistent with zero: Act = —0.27" ± 0.32" 
and AS = 0.14"±0.39". Finally, the median ratio of fluxes 
of bright objects detected in our radio map to the fluxes 
reported in the NRAO VLA Sky Survey (NVSS; Condon 
et al. 1998) is - 1.19±g;le, confirming the accuracy of the 
flux calibration. 

However, our detection is inconsistent with a non- 
detection at 2.1 GHz presented by Hatsukade et al. (2012, 
S2.1 < 34.2 /xJy at 3cr) as it implies an extremely steep 
(and unphysical) spectral index < —5.9. Extrapo- 

lating from our 1.43 GHz detection, the expected signal at 



2.1 GHz would be ~ 250 ^Jy assuming a = —0.75. How- 
ever the equatorial declination of GRB 021211 of +6°44' 
makes it difficult to observe it using east-west arrays, such 
as ATCA, as the beam is highly elongated, i.e. 2" x 51" at 
2.1 GHz. This may pose some problems in the detection 
of sources, and indeed in our VLA map we have found 
two additional sources 23 with 1.43 GHz fluxes of ~ 650 
and 350 ^iJy, respectively, which are not detected in the 
2.1 GHz ATCA map of Hatsukade et al. (2012). Further 
observations at various radio frequencies are needed to re- 
solve this issue. 

Our 1.43 GHz detection is consistent with the 
(sub)millimctcr limits of Smith et al. (2005) and Prid- 
dey et al. (2006). Namely, they did not detect the host 
of GRB 021211 at 850 /im (0.3 ± 1.9 mJy) and 1.2 mm 



With the following radio R.A. and deck: 08:09:01.315, +06:43:03.91; 08:09:12.842, +06:43:54.53; and 08:09:11.746, +06:44:05.87. 
With the following VLA R.A. and deck: 08:09:10.992, +06:41:26.20; and 08:09:31.704, +06:41:58.00. 
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(0.07 ± 0.53 mJy), respectively. The 850 /im limit implies 
a submillimcter-to-radio spectral index a 359 < 0.53, con- 
sistent with most of the models presented by Carilli & Yun 
(1999, 2000, their Figures 1 and 3, respectively). Assuming 
spectral energy distribution (SED) templates of Arp 220 
and M82 (Silva et al. 1998), the 1.2 mm limit corresponds 
to a 3cr limit of SFR < 700-1100 A^yr^ 1 (for the 850 
limit these estimates are ~ 30% higher), consistent with 
our radio SFR ~ 825 M yr" 1 (Table 2). However, this 1.2 
mm flux limit corresponds to an SFR ~ 100-500 M yr" 1 
for many other SED templates (Silva et al. 1998; Iglesias- 
Paramo et al. 2007; Michalowski et al. 2008, 2010a). 
Hence, deeper (sub)millimctcr observations (rms ~ 0.1- 
0.2 mJy) are needed to verify if our radio detection is 
inconsistent with (sub) millimeter data, which would in- 
dicate a significant AGN contribution to the radio flux of 
the host of GRB 021211. 

5.2. Star Formation Rates of the GRB Host Population 

The SFRs of GRB hosts are shown as a cumulative dis- 
tribution on Figure 3. The high-SFR boundaries were 
calculated using the radio detections and upper limits, 
whereas the low-SFR boundaries were obtained by sub- 
stituting the radio SFR upper limits with the lower limits 
from the UV. We found that at least - 63% (M5/24) 24 
of all our GRB hosts at z < 1 have SFR < lOOMgyr" 1 
and only < 8% (<2/24) 25 could have SFR > 500 A^yr" 1 . 
This implies that it is rare (< 33% chance, <8/24) 26 for 
a GRB to reside in an ULIRG. This is consistent with 
the contribution of ULIRGs to the cosmic star formation 
history being < 10% at z < 1 (Lc Floc'h et al. 2005). 

Even though high star-forming GRB hosts are rare at 
z < 1, the SFR of GRB 021211 alone constitutes as much 
as - 22% of the summed SFR of all z < 1 GRB hosts, 
even when we sum over radio upper limits (and hence its 
contribution is higher in reality). Hence, such high star- 
forming GRB hosts likely dominate the contribution of 
this population to the cosmic star formation history. 

The average radio SFR of GRB hosts can be assessed 
using the average radio flux of the GRB hosts undetected 
in our radio observations. For each host we converted the 
flux at the GRB position at the observed frequency to that 
at the rest-frame 1.43 GHz, using a radio spectral index 
of —0.75. In this way we obtained a weighted mean of the 
flux equal to — 13 ± 16 /iJy. Hence, we did not detect the 
GRB host population even when averaging the data. At 
least such level of rms has to be reached in future GRB 
host surveys to obtain significant number of detections. 
At the mean redshifts of these hosts, z = 0.53, this cor- 
responds to a 3er upper limit of SFR ra di < 15Af yr _1 . 
Hence, the general population of GRB hosts is below the 
LIRG limit (L m < lO n L , or SFR < 17.2 A^yr" 1 us- 
ing the conversion of Kennicutt 1998). It is expected that 
LIRGs do not dominate our GRB host sample, because 
LIRGs dominate the cosmic star formation history only 
above z ~ 0.7 (with their contribution rising to ~ 70% at 

24 > 50% (>5/10) for TOUGH subset only. 

25 0% (0/10) for TOUGH subset only. 

26 < 50% (<5/10) for TOUGH subset only. 

27 Assuming fourth octile of water vapor; 

http : //almascience . eso . org/call-f or-proposals/sensitivity-ca 

28 http : //tao . it . swin. edu. au/mock-galaxy-f actory/ 



z = 1, Le Floc'h et al. 2005; and staying at this level at 
least up to z ~ 2.3, Magnclli et al. 2011). 

The full ALMA with 50 antennas will reach an rms 
sensitivity of ~ 0.023 mJy at 345 GHz in 1 hr 27 . This 
corresponds to SFR ~ 5-20M Q yr _1 at z = 10 (using 
SED templates of Silva et al. 1998; Iglesias-Paramo et al. 
2007; Michalowski et al. 2008, 2010a), so ALMA will eas- 
ily detect GRB hosts basically at any redshift within a 
few hours, because the UV lower limits on SFRs are of the 
order of - 1 A/ yr" 1 (Table 2). 

To summarize, the overall picture is that z < 1 GRB 
hosts have modest SFRs (as suggested by Stanway et al. 
2010), but a small fraction (~ 4-8%) of them have un- 
dergone an extreme star formation episode. However the 
latter claim suffers from poor number statistics. 

5.3. The Relation to Other Galaxies: Do GRBs Trace 
Star Formation in an Unbiased Way? 

In order to investigate whether the GRB host population 
is consistent with the general population of star-forming 
galaxies at similar redshifts, we show their SFR distribu- 
tion on Figure 3. The comparison to other galaxies must 
be done carefully, because the probability that a galaxy 
with given SFR is included in a usual galaxy sample de- 
pends only on the number density of such objects (as long 
as this SFR corresponds to a flux higher than the sample 
selection threshold). This is not the case for a GRB host 
sample, because, assuming that GRBs trace star forma- 
tion in an unbiased way, a galaxy with higher SFR is more 
likely to host a GRB and, in turn, to be selected into the 
GRB host sample (e.g. Natarajan et al. 1997; Fynbo et al. 
2001). In order to account for this, we weighted the cumu- 
lative distributions of other galaxies by their SFRs, i.e. the 
curves for other galaxies correspond to the fraction of total 
star formation in the sample contributed by galaxies with 
SFRs lower than a given SFR. 

It is apparent from Figure 3 that the SFR distributions 
of GRB hosts and of simulated galaxies at z = 0.51 (Cro- 
ton et al. 2006) 28 , produced in a semi-analytical model and 
based on the Millenium simulation (Springel et al. 2005), 
are fully consistent. 

Similarly, the SFR distribution of z < 1 GRB hosts is 
consistent with that of z ~ 1 Lyman break galaxies (LBGs; 
from Basu-Zych et al. 2011, SFRs from 24 /xm and rest- 
frame UV photometry) and z ~ 0.84 Ha emitters (HAEs; 
from Sobral et al. 2009, SFRs from 24 /mi photometry and 
Ha fluxes). We note that the median stellar masses of 
these LBGs (M* ~ 10 9 5 Af o ; Basu-Zych ct al. 2011) and 
HAEs (M* ~ 1O 1O1 M ; Sobral et al. 2011) are also con- 
sistent with that of GRB hosts (M* ~ lO 9 - 3 " 9 - 7 M ; Cas- 
tro Ceron et al. 2010; Savaglio et al. 2009). Moreover, as 
shown in Section 5.4, the dust attenuation we derived for 
GRB hosts is consistent with that of LBGs and HAEs. 

An apparent inconsistency at low SFRs of the GRB host 
samples with the LBG and HAE populations (the GRB 
population extends to lower SFRs) is an effect of higher 
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Fig. 3. — Cumulative distribution of SFRs of GRB hosts in the z < 1 TOUGH unbiased (red area) and the GRB-SN (blue area) subsets. The 
high-SFR boundaries (colored solid lines) are constructed using the detections and limits of SFR ra( jj (Table 2), whereas the low-SFR bound- 
aries (colored dotted lines) are constructed using the SFRtjv for galaxies not detected in the radio. The upper rr-axis gives the corresponding 
infrared luminosity according to SFR(M yr" 1 ) = 1.72 X 1O~ 1o Ltr(L ) (Kennicutt 1998). Wc found that at least ~ 63% (>15/24) of GRB 
hosts at z < 1 have SFR < 100 Mq yr _1 and only < 8% (<2/24) could have SFR > 500 Mq yr" 1 . For comparison, the SFR distributions of 
z = 0.51 simulated galaxies (Croton et al. 2006), z ~ 1 Lyman break galaxies (LBG; Basu-Zych et al. 2011), z ~ 0.84 Ha emitters (HAE; 
Sobral et al. 2009), and z <~ 2-3 submillimeter galaxies (SMGs; Michalowski et al. 2010a, b) are shown (labelled lines, of which the right lines 
represent dust-corrected SFRs). These distributions were weighted by SFR (so they reflect the fraction of total star formation in the sample 
contributed by galaxies with SFRs lower than a given SFR) to allow a comparison with the GRB host population, which is likely selected 
based on SFRs (see Section 5.3). It is evident that current SFR limits imply that the GRB host population is consistent with star-forming 
galaxies at similar redshifts (simulated, LBGs, and HAEs) and is inconsistent with SMGs. Since we did not detect most of the targets, the 
distributions of the z < 1 TOUGH unbiased and GRB-SN subsets are consistent (the overlap of the blue and red areas is significant). 



flux detection threshold for the latter. Namely, the lim- 
iting magnitude for LBGs of u < 24.5 mag corresponds 
to SFR > 0.5A/ Q yr _1 , whereas the limiting luminosity 
for HAE of L Hq > 10 415 erg s -1 corresponds to SFR > 
2.5 M Q yr -1 . Hence, galaxies with SFR < 0.5 M Q yr -1 are 
not present in the LBG and HAE samples, because they 
are below the detection limits. Indeed when we restricted 
the simulated galaxies (which are consistent with the GRB 
host population) to galaxies above these limits, their dis- 
tributions are consistent with those of LBGs and HAEs. 

A Kolmogorov-Smirnov (K-S) test resulted in a proba- 
bility of - 15% that the UV SFRs of all our GRB hosts 
and z < 1 LBGs are drawn from the same population. 
However, for z ~ 0.84 HAEs such probability is negligible 
(~ 10~ u ), showing that HAEs have systematically higher 



UV SFRs than GRB hosts. Figures 2 and 3 show that 
our current limits on the radio SFRs of GRB hosts are not 
deep enough to test whether the total SFRs of HAEs are 
also higher than that of GRB hosts. 

As shown in Figure 3, the SFRs of the z < 1 GRB 
hosts are clearly inconsistent with those of submillimeter 
galaxies (SMGs), dusty high star- forming z ~ 2-3 galaxies 
(Michalowski et al. 2010a, b, SFRs from total infrared emis- 
sion and rest-frame UV photometry), even when only radio 
upper limits for GRB hosts arc taken into account (colored 
solid lines). This is also expected from the fact that GRB 
hosts are much less massive (M* ~ 1O 9 3_9 ' 7 M0; Cas- 
tro Ceron et al. 2010; Savaglio et al. 2009) than SMGs 
(M* ~ ioio-4-ii.3 Mq . Borys et al 2005; Michalowski 
et al. 2010a, 2012; Hainline et al. 2011; Bussmann et al. 
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2012; Yun ct al. 2012). We note however, that the stellar 
mass estimates for GRB hosts are based on samples biased 
against dusty galaxies, so the general population of GRB 
hosts may include galaxies as massive as SMGs. 

The inconsistency of the z < 1 GRB host and SMG pop- 
ulations is also revealed by the K-S test giving a negligible 
probability of ~ 10~ 8 that the UV SFRs of GRB hosts and 
SMGs are drawn from the same population. If we restrict 
the analysis to z < 1 SMGs then the probability increases 
to ~ 63% (consistent samples), but these z < 1 SMGs 
have median SFRir ~ 60Mgyr , very close to our lim- 
its of SFR ra dio for GRB hosts, so if GRB hosts are similar 
to z < 1 SMGs, then the majority of them would need to 
be just below the detection limits, which is unlikely 

We note that the comparison between GRB hosts, 
SMGs, LBGs, and HAEs involves SFRs derived from the 
radio, 24 //m, total infrared, UV, and Ha luminosities, but 
we do not expect any systematic offset between these esti- 
mates. Namely, Elbaz et al. (2010, their Figure 2) showed 
that the observed 24 /im luminosity is correlated with the 
total infrared luminosity, which, in turn, is correlated with 
the radio luminosity (Condon 1992). Similarly, Wijcsinghc 
et al. (2011, their Figure 8) showed that SFRs from Ha 
and UV are well correlated. 

To summarize, our data allow a significant range of 
SFRs for z < 1 GRB hosts, so their distribution is con- 
sistent with that of the general population of galaxies (at 
least with LBGs and simulated galaxies) having a median 
SFR of a few Mgyr . If this conclusion is confirmed in 
deeper radio observations, it would indicate that GRBs 
trace a large fraction of the overall star formation, and are 
therefore less biased indicators than once thought. On the 
other hand, if deeper radio observations reveal that the 
total SFRs of GRB hosts are very close to their UV SFRs, 
then GRB hosts will not be consistent with tracing the 
overall population of star-forming galaxies. 

Indeed, there are some indications that GRB hosts are 
different than other galaxies at similar redshifts. Previ- 
ous studies have showed that GRB hosts are metal-poor 
(Fynbo et al. 2003; Prochaska et al. 2004; Gorosabel et al. 
2005; Sollerman et al. 2005; Hammer et al. 2006; Stanek 
et al. 2006; Wiersema et al. 2007; Christensen et al. 2008; 
Modjaz et al. 2008, 2011; Levesque et al. 2010a, b; Leloudas 
et al. 2011), fainter, and more compact than SN hosts, 
as well as that GRBs themselves are much more concen- 
trated in the UV-bright regions of their hosts than SNe 
(Fruchter et al. 2006; Bloom et al. 2002). Moreover, Wain- 
wright ct al. (2007, their Figures 3 and 6) found that in 
the z < 1 GRB host population the ratio of irregular to 
regular galaxies is 2:1, different than 1:2 for field galax- 
ies; and that GRB hosts are a factor of two smaller than 
field galaxies. This suggests a potential bias in the GRB 
host population (though all these studies were based on 
optically-biased GRB host samples). 

However, Leloudas et al. (2010, their Figure 3) showed 
that the distribution of GRBs within their hosts is in fact 
consistent with that of type-Ic SNe and of Wolf-Rayet 
stars, whereas Svensson et al. (2010, their Table 3 and 
Figure 8) did not find any significant difference between 
the absolute magnitudes of GRB and SN hosts. More- 
over, Fynbo et al. (2006, their Figure 2) and Savaglio 

29 > 80% (>8/10) for TOUGH subset only. 



et al. (2009, their Figure 18) showed that the metalic- 
ity of GRB hosts is consistent with (or even higher than; 
see also Savaglio et al. 2012) that of damped Lya sys- 
tems at corresponding redshifts. Similarly, Levesque et al. 
(2010c) found that some GRBs explode in metal-rich envi- 
ronments. Moreover, optically-dark GRBs, missed in most 
previous studies, were claimed to be hosted in more dusty 
and metal-rich hosts (Perley et al. 2009; Greiner et al. 
2011; Hunt et al. 2011; Kriihler et al. 2011; Svensson et al. 
2012). Finally, Conselice et al. (2005) found that GRB 
hosts are smaller than field galaxies, but only at z < 1.2; 
and that these two populations at z < 1 are consistent with 
regards to their concentrations and asymmetries. Hence, 
the biases of the GRB host sample in terms of morphology 
and metallicity are far from being well understood. 

Regarding observational biases, Fynbo ct al. (2009) 
showed that the GRB host sample selected based on a 
requirement of an optical spectroscopic redshift is not rep- 
resentative for all GRBs, likely biased against dusty ob- 
jects. Moreover, SNe exploding in (U)LIRGs are expected 
to be highly extinguished by dust (Mattila et al. 2004, 
2007, 2012; Melinder et al. 2012), implying that samples 
of optically-selected GRBs may miss completely those ex- 
ploding in high star-forming galaxies. However, Figure 3 
does not reveal any differences between the radio prop- 
erties of the optically biased pie-Swift sample of GRB- 
SN hosts and those of the z < 1 TOUGH unbiased sub- 
set. Namely, the high-SFR bounds of the distributions 
are similar. This is mainly because we did not detect the 
majority of the targets and the limiting depths were sim- 
ilar. Similarly, the K-S test resulted in 40% probability 
that the SFR UV for the TOUGH and GRB-SN subsets 
are drawn from the same parent population. Larger sam- 
ples and deeper radio data are necessary to investigate this 
issue. 

5.4. Ultraviolet Attenuation 

The cumulative distribution of the dust attenuation is 
shown on Figure 4. It should be seen as an upper limit, 
because it is based mostly on radio non-detections. We 
found that > 88% (>21/24) 29 GRB hosts at z < 1 have 
Atjv < 6.7 mag, i.e. Ay < 3 mag assuming a Small Mag- 
ellanic Cloud (SMC) extinction curve (Atjv = 2.22 Ay; 
Gordon et al. 2003). This is consistent with the previous 
observational (Kann et al. 2006; Schady et al. 2007, 2010, 
2012; Kann et al. 2010; Savaglio et al. 2009; Han et al. 
2010; Liang & Li 2010; Greiner et al. 2011; Watson 2011; 
Watson & Jakobsson 2012; Zafar et al. 2010, 2011) and 
theoretical (Lapi et al. 2008; Mao 2010) results that GRB 
hosts are weakly obscured by dust with very few excep- 
tions (Tanvir et al. 2004, 2008; Castro Ceron et al. 2006; Le 
Floc'h et al. 2006; Michalowski et al. 2008; Kiipcii Yoldas, 
et al. 2010; Hunt et al. 2011; Kriihler et al. 2011; Svensson 
et al. 2012). This result is also confirmed for the entire 
TOUGH sample (Hjorth et al. 2012), which is not biased 
against dusty systems. Our finding is also consistent with 
a fraction of 5-36% of core collapse SNe in normal galaxies 
with A v > 3.7 mag (Section 2.4 of Mattila et al. 2012). 

It has been claimed that dust is responsible for the 
optical faintness of the so-called dark GRBs, i.e. those 
with optical-to-X-ray spectral index /3qx < 0.5 (Ramirez- 
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Fig. 4. — Cumulative distribution of ultraviolet (and optical) dust attenuation of GRB hosts in the z < 1 TOUGH unbiased (red solid 
line) and GRB-SN (blue solid line) subsets. Since the distributions arc mostly based on radio non-detections they should be taken as upper 
limits to Atjv, i.e. the area to the right of the lines is ruled out by the data. The upper x-axis was derived by assuming an SMC extinction 
curve (A uv = 2.22 Ay; Gordon et al. 2003). We found that > 88% (>21/24) of the GRB hosts have A uv < 6.7 mag, i.e. A v < 3 mag. For 
comparison, the distributions of z ~ 1 Lyman break galaxies (LBG; Basu-Zych ct al. 2011), z ~ 0.84 Ha emitters (HAE; Villar et al. 2008; 
Sobral et al. 2009), and z ~ 2-3 submillimcter galaxies (SMGs; Michalowski et al. 2010a, b) are shown (labelled black lines). The attenuation 
at Ha was converted to that in the ultraviolet by assuming an SMC extinction curve (Atjv = 1-78 An a ; Gordon et al. 2003). The distributions 
for GRB hosts are consistent with those of LBGs and HAEs, yet likely inconsistent with that of SMGs. 



Ruiz et al. 2002; Jakobsson ct al. 2004; Perley et al. 2009; 
Greiner et al. 2011; Kriihler et al. 2011). None of our z < 1 
GRBs were dark (Jakobsson et al. 2004; Butler et al. 2005; 
Fynbo et al. 2009; Xu et al. 2009, though there is no X-ray 
data for four pre-Swift GRBs in our sample). Hence, none 
of the afterglows of the GRBs in our z < 1 sample seems 
to be particularly dust-obscured (though the effect of dust 
at z < 1 is less severe than at higher redshifts, i.e. for 
the same amount of dust a GRB at higher redshift may 
be classified as "dark", because the observed-frame opti- 
cal emission corresponds to the UV wavelengths strongly 
affected by dust). This is consistent with low levels of 
dust attenuation we find for their hosts, but we note that, 
in principle, the spatially-integrated attenuation of a host 
reported here may be inconsistent with the line-of-sight 
extinction derived from an afterglow, as the latter probes 
dust distributed only over a narrow opening angle. 

GRB 021211, with a radio-detected host, was initially 
called "dark" due to its optical faintness (Crew et al. 2003; 



Fox et al. 2003; Li et al. 2003), but the lack of X-ray data 
makes it impossible to classify it according to the more 
rigorous /?ox definition. However, Fox et al. (2003) found 
that the optical afterglow was not severely reddened, so 
the GRB must have occurred away from the highly dust- 
attenuated regions suggested by our radio detection. 

We also note that only < 50% of all our GRB hosts at 
z < 1 could have Ay > 2 mag (Atjv > 4.4 mag), i.e. sim- 
ilar dust attenuation levels to that in SMGs (Smail et al. 
2004; Swinbank et al. 2004; Borys et al. 2005; Michalowski 
et al. 2010a,b; Hainline et al. 2011). Indeed, in Fig- 
ure 4 we show that the Atjv distribution for SMGs (from 
Michalowski et al. 2010a, b) displays high attenuation lev- 
els, very likely inconsistent with that of z < 1 GRB hosts, 
given that the lines for GRB hosts represent upper lim- 
its on dust attenuation. This is consistent with a ten- 
dency that the presence of a GRB typically selects dwarf 
galaxies that are generally less dusty. However, the most 
star-forming GRB hosts may contain significant amounts 
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of dust, comparable to those of SMGs (as suggested by 
Michalowski et al. 2008). 

In Figure 4 we also show the attenuation distribution 
of z ~ 1 LBGs from Basu-Zych et al. (2011). Our GRB 
host distribution is mostly based on upper limits on A\jy, 
but it is clear that these two distribution may be consis- 
tent, i.e. the fraction of objects with very low dust atten- 
uation is similar (~ 20-30%) and the LBG distribution 
is always below that of the GRB hosts. Deeper radio 
or far-IR data would be necessary to confirm that these 
samples are truly consistent. If this is the case, then we 
can expect a median Atjv ~ 1-6 mag for GRB hosts and 
therefore that their SFR ra( jio should be, on average, a fac- 
tor of ~ 4 higher than their SFRtjv- Assuming a typical 
SFR UV = lMoyr" 1 , then the SFR radio =4M yr" 1 cor- 
responds to the observed- frame 1.43 GHz fluxes of 33, 8, 
and 3 /iJy at z = 0.3, 0.6, and 1.0, respectively. Hence, if 
the attenuation of GRB hosts is indeed similar to that 
in LBGs, then current radio interferometers (including 
EVLA) will struggle to detect them beyond z = 0.6. 

For HAE we converted the attenuation at Ha to that 
at the UV assuming an SMC extinction curve (j4uv = 
1.79 Aft a ; Gordon et al. 2003). The comparison with the 
dust attenuation of GRB hosts and HAEs is less clear. The 
sample of Villar et al. (2008) contains much fewer galaxies 
with very low attenuation compared to the z < 1 GRB 
host sample, but this is not the case for the sample of So- 
bral et al. (2009). In any case, in order for a GRB host 
sample to be consistent with the HAE sample, the attenu- 
ation of the hosts with current limits at A\jy ~ 5 — 7 mag 
would need to be very close to these limits. 

5.5. Radio Spectral Indices 

The radio spectral index of the host of GRB 980425 
turned out to be very steep (af'| 4 < —1.44; Michalowski 
et al. 2009) and was interpreted as a sign of the domi- 
nant old stellar population. Similarly, for the GRB 021211 
host, using our detection at 1.43 GHz and the upper limit 
of 45 jttJy at 8.46 GHz reported by Fox et al. (2003) we 
derive a steep af;|| < —1.12. 

On the other hand, the spectral index of the GRB 
031203 host is flatter a\ f 9 = -0.53 ± 0.50. This suggests 
a higher contribution of free-free emission (or synchrotron 
self-absorption; Condon 1992) and, hence, younger stellar 
population (Brcssan et al. 2002; Cannon & Skillman 2004; 
Hirashita & Hunt 2006; Clemens et al. 2008). This value 
is also consistent within errors with the spectral slopes of 
a ~ —0.7 to —0.8 found for star-forming galaxies both lo- 
cal and at high redshifts (Condon 1992; Dunne et al. 2009; 
Ibar et al. 2010). 

6. CONCLUSIONS 

We present radio continuum data for a sample of 30 
GRB hosts including 22 new observations. We detected 
three targets. The derived limits on the SFRs show that 
at least ~ 63% of the GRB hosts have SFR < 100 M yr" 1 



and that at most ~ 88% of GRB hosts have Atjv < 6.7 
mag, i.e. Ay < 3 mag. The average flux of non-detected 
hosts at z ~ 0.5 sets an upper limit of SFR < 15 Mq yr _1 . 
Using our radio data in conjunction with the rest-frame ul- 
traviolet data we found that the distributions of SFRs and 
ultraviolet attenuations of GRB hosts are consistent with 
those of other star-forming galaxies at z < 1. This is con- 
sistent with the hypothesis that GRBs trace cosmic star 
formation, but further studies of morphology and metal- 
licitics of GRB hosts are required to understand potential 
biases in this sample. 
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